Quantum Monte Carlo simulations of the t-Jz model with stripes on the square lattice 
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Finite temperature quantum Monte Carlo simulations are performed on the anisotropic t — J 
model and in particular on its Ising limit. Straight site-centered stripes are imposed by an on- 
site potential representing external mechanisms of stripe formation. In this model, we show that, 
even though charge inhomogeneity exists at a high temperature, the anti-phase ordering of the spin 
domains between stripes occurs at a much lower temperature. The magnetic correlations at this 
spin ordering crossover are analyzed. The stripes show metallicity, with absence of hole attraction. 
Comparison between this model and others that have been proposed to explain or describe stripes, 
as well as possible relations with experimental features on underdoped cuprates are discussed. 
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The nature of the underdoped region and of the pair- 
ing mechanism in high-T c cuprates is still a matter of 
strong controversies. From the experimental point of 
view, the relationship between the pseudogap (PG) and 
the superconducting gap (SG) is still under intense dis- 
cussion. Older data suggested a pseudogap smoothly 
connecting with the. .superconducting gapu while some 
more recent studiestm emphasize a different origin and 
behavior of the PG and the SG, the PG ending in- 
side the superconducting phase at a "quantum critical 
point" D From the theoretical point of view, the stripe 
scenario, which is an almost unavoidable consequence 
of neutron experiments results on La2-xSr x Cu04cllj and 
perhaps on YE^CusOg+sB, offers a natural explanation 
for PG. Consistently with the experimental controversy 
about the nature of PG, there are widely diverging views 
about the origin of stripes and its relationship with su- 
perconductivity. In the theory by Emery, Kivelson and 
coworkers!], the pseudogap is identified with the SG and 
stripes are a key ingredient to explain superconductivity 
in the cuprates. In another view,LI stripes are regarded 
as competing with a uniform gas of hole pairs and hence 
with superconductivity. In szet another approaches, the 
stripes exclude hole pairing]!!] There are even more im- 
portant differences regarding the origin of stripes. In 
Rcf. ||, charge inhomogeneity appears as a combined ef- 
fect of phase separation and long-range Coulomb repul- 
sion and the stripe phase can be thought pas a Wigner 
crystal. According to White and ScalapinoO, stripes are 
already present in the simple t — J model at physical val- 
ues of the parameters. In Ref. [l0| the driving mechanism 
is the formation of strong singlets across a hole. In this, 
sense, the stripes may be regarded as domain walls.Ell 
On the other hand, in some other views, the stripes are 
not inherent to two-dimensional (2D) extended t — J 
models buJ_are due to, for example, electron-phonon 
couplingEMij Following these alternative views, in this 
paper we formulate a model assuming that the origin of 
stripes in the cuprates is non-intrinsic to the 2D elec- 
tronic correlations and our main goal is to examine the 



physically relevant properties of such model. 

It is also well known the difficulty of studying micro- 
scopic models of correlated electrons either by analytical 
or numerical techniques. The already mentioned work 
on variants of the t — J model by using the DMRG 
methocoij is affected by limitations of this method, par- 
ticularly the open boundary conditions adopted. The 
widely used Lanczos diagonalization could not deal with 
clusters large enough to reproduce the charge inhomo- 
geneities. In the present study, we use the conventional 
finite temperature quaptum Monte Carlo (QMC) method 
(world-line algorithm)!!^ which allows the study of rea- 
sonable large clusters with fully periodic boundary condi- 
tions. As it is well-known, QMC simulations of fermionic. 
models are affected by strong "minus sign problem"E£l 
that virtually makes impossible these kind of studies at 
very low temperatures. However, as we discuss below, in 
the model we consider the "minus sign problem" is not 
as much severe as in the plain 2D t — J model, and hence 
we are able to look at physical features which appear at 
not too low temperatures. 

Let us first introduce the model here studied, l-^tue 
first place, as in many other studies on this subjectE3~o, 
we impose the presence of straight site-centered stripes 
by an on-site potential. This on-site potential repre- 
sents the effects of Coulomb potential due to out-of-plane 
ions, electron-phonon coupling, a-b plane anisotropy or 
other structural details, etc. Since in-plane long-range 
Coulomb repulsion can not be included in QMC simu- 
lations, its effects can eventually also be represented by 
the on-site potential. The confinement of holes to the 
stripes strongly reduces the "minus sign problem" . As 
we show below, this problem strongly constrains the val- 
ues of the on-site potential accessible in the simulations. 
A further alleviation of this problem comes from reduc- 
ing quantum spin fluctuations. We are thus lead to an 
anisotropic t — J modelr-and in the Ising limit to the 
so-called t— Jz model£iH23 Again, the "minus sign prob- 
lem" strongly constrains the off-diagonal exchange term 
we can deal with. On the other hand, this Ising limit 
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is not only convenient from the numerical point of view 
but has an additional value, specially if one is interested 
in finite temperature effects. In fact, a three-dimensional 
(3D) AF (short- or long-range) order would imply, at a 
mean-field level, a staggered field acting on the 2D t-J 
model and in turn would induce an enhancement of the 
z^-component of the Heisenberg term as a second order 
process. 

The Hamiltonian of the anisotropic t — J model is: 



Hu = - 1 Y, (4^ + h.c.) 



jJ2(~(StSr + SrSf) + S*S* - -n^) (1) 



(ij) 

where the notation is standard. , 
by an effective on-site potential:! 



,Jhe stripes are induced 



H s tr — 



E 



(2) 



e si = — 2e s < (2e s ) for sites on (outside) the stripe. 
Then, the total Hamiltonian of our model is H = H t j + 
H str - Most of our calculations were performed on an 8 x 8 
cluster with fully periodic boundary conditions and with 
eight holes which corresponds to a filling of x = 1/8. The 
imposed stripes involves columns separated by three-leg 
ladders, as in the original picture in Ref. [5| 

Our simulations were done at J/t — 0.35, a value 
generally accepted to describe the cuprates, and also at 
J/t = 0.7. In this second case, we have seen essentially 
the same physical behavior but at a higher temperature 
scale and with a milder minus sign problem. As usual, t 
is chosen as the unit scale of energy and temperature. e s 
was varied between 0.3 and 2.0, and 7 = 0.0 (the Ising 
limit), 0.25 and 0.5 were examined. In the limit e s —> 00 
each stripe would be at quarter filling. Simulations were 
also performed for the 12 x 12 cluster with twelve holes 
and also two equidistant stripes. This corresponds to a 
smaller density x = 1/12 and the spin domains in be- 
tween the stripes are five-leg ladders. 

The QMC algorithm employed is a straightforward ex- 
tension of the world-line cme successfully used to study 
the 2D Heisenberg modeLEJ Besides the cube and pla- 
quette local moves, we have kept the global moves that 
change the total S z of the system. Most of the calcu- 
lations were performed with r = (3/M = 0.083, where 
(3 = 1/T and M is the Trotter number. The average of 
the sign of exp — j3H, is shown in Fig. |l](a). It can be 
seen that {sgn} is smaller for larger 7 and for smaller 
e s . In few words, the more isotropic in the spin space 
and the more homogeneous is the hole movement in real 
space, the worse becomes the minus sign problem. In 
Fig. [y(b) , the hole occupancy of sites on the stripes are 
shown as a function of temperature for various sets of 
parameters (J,e s , and 7). As expected, the stripe filling 
increases monotonically as T decreases. At large temper- 
ature there is already a filling of the stripes larger than 



the nominal one (x = 1/8) and it is driven solely by the 
on-site potential. As the temperature goes to zero, the 
stripe density saturates at a value smaller than 1/2 and 
it is apparent a small dependence on J, i.e., there is a 
correction due to the magnetic correlations of the t — J 
model. In Fig. 0(c), the hole density profile is shown for 
J = 0.35 and several values of e s at the lowest temper- 
atures reached. The sharpness oi-these profiles could be 
measured by neutron scattering.Ea 
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FIG. 1. (a) Averaged sign in QMC simulations on the 8x8 
cluster with eight holes. The sets (J,e s , 7) are indicated on 
the plot. Stars correspond to 12x 12, twelve holes (x = 0.083). 
(b) Hole density on the stripes vs. temperature, 8x8 cluster. 
The dotted line corresponds to the uniform hole distribution 
(x = 1/8). (c) Hole density profile (I — and 4 correspond to 
the stripes) for the t — Jz model with J = 0.35, T < 0.1 and 
for several values of e s as indicated on the plot (T = 0.15 for 
e s = 0.3). 

We start now to show the main features observed.il 
Results of the computation of magnetic and charge static 
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structure factors, S(k) and C(k) (Fourier transformed 
spin-spin and hole-hole correlations functions respec- 
tively) are partially summarized in Fig. ^(a). In the 8x8 
cluster (x — 1 /8) , at a rather large T there is a crossover 
in the peak of the structure factor from (k x , k y ) = (tt, 0) 
(with the x-axis perpendicular to the stripes direction) 
to (25,0) with 5 — 7r/4. This crossover, together with 
the behavior shown in Fig. |l|(b) indicates a non-trivial 
behavior of the charge ordering. One is tempted to_t£rm 
this crossover as the "charge ordering" temperaturej^j al- 
though this concept is somewhat arbitrary in our model. 
The most important feature is that at a temperature 
much lower that this crossover there is a second crossover 
in the spin sector signalled by a change in the peak of the 
magnetic structure factor from (tt, tt) to (tt ± S, tt). This 
peak very much resembles the one observed in neutron 
scattering experimentaj'Ll signalling the presence of the 
"incommensurate phase" in underdoped cuprates. Fol- 
lowing Ref. ||, we call T 2 this lower crossover at which 
a "spin ordering" occurs. A similar behavior is observed 
in the 12 x 12 lattice with 12 holes (x = 0.083). Follow- 
ing the behavior observed in underdoped cuprates, T 2 * 
occurs at a temperature higher than the one observed 
for similar parameters but for x — 0.125. The peaks of 
S'(k) and C(k) are the same as above except that for this 
smaller hole doping, 5 — tt/6. A more detailed evolution 
of the peaks with temperature is shown in Fig. ^](b). In 
particular, the crossover in the charge sector is rather 
smooth. In the spin sector, the weight of the (tt, tt) peak 
is strongly reduced below the crossover, but S(tt ± 5, tt) 
is definitely nonzero above it. 
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FIG. 2. (a) Phase diagram in the temperature-on-site 
potential plane. The curves at high (low) T correspond to 
a crossover in the charge (magnetic) structure factor. The 
sets (J, 7) are indicated on the plot. The star corresponds 
to 12 x 12, x = 0.083, J = 0.7, 7 = 0.0. (b) Magnetic and 
charge static structure factors vs. temperature for J/t — 0.35, 
[e s | =0.5, 7 = 



We have detected important hysteresis effects in the 
crossover region. The results shown in Fig. |2|(b) were 
obtained by starting the simulation from a configuration 
with 7r-shifted spin domains. If the starting configuration 
consists of in-phase spin domains the crossover tempera- 
ture is pushed to smaller values. In this case, for many 
values of the parameters J, e s and 7 (specially for large 
J and e s ) we could not detect that crossover to the low- 
est temperature attainable. In most cases, starting from 
a randomized spin configuration, behaviors like those of 
Fig. [2] are recovered but we cannot rule out the possibility 
of in-phase domains being rnacc stable than or degener- 
ate with anti-phase domainsQEJ at low temperatures and 
for very large J and e s . 

We can analyze the crossover in the spin sector at 
a more microscopic level by looking at the real-space 
spin-spin correlations that experience the most impor- 
tant changes at this crossover. These correlations are 
between sites two lattice spacings apart in the same row 
across a stripe (labelled Si), and between sites belonging 
to the center leg of two consecutive three-leg ladders in 
the same row (S2). For completeness we have also com- 
puted the correlations between nearest neighbor (NN) 
sites in the center leg of a three- leg ladder (S3), and be- 
tween sites at the maximum distance along this leg (S4). 
S5 and Sg are the correlations between NN and next NN 
(NNN) sites along the stripes. The spin-spin correlations 
have been normalized in such a way that their maximum 
(minimum) value is +1 (—1) for the z-components of the 
two spins fully aligned or ferromagnetic (FM) (respec- 
tively anti-aligned or AF). In Fig. ||(a), corresponding 
to J/t = 0.35, e s = 1.0 it can be seen that the corre- 
lations Si and S2 across the stripe are positive at high 
temperature and they increase as T is lowered. Around 
T « 0.12 these correlations reach their maximum value 
and as the temperature is further lowered they suddenly 
become fully AF and remain negative down to the low- 
est temperature reached. These changes from FM to AF 
indicate that the intervening mostly hole-free ladders are 
in-phase (anti-phase) above (below) T w 0.11, S3 and 
Si (only shown in Fig. ^(a)) show a monotonic behav- 
ior as the temperature is decreased. They indicate a full 
polarization of the spin domains at low temperatures. A 
similar behavior can be observed for smaller e s and larger 
J (Fig. |(b)). For the e s = 1.0 but J = 0.7 (Fig. |(c)), 
T 2 * increases and it increases further for smaller density 
(x = 0.083 on the 12 x 12 cluster). The same features 
survive when a XY term in the exchange interaction is in- 
cluded, as it can be seen in Fig. ||(d) for 7 = 0.25 and 0.5. 
As the Heisenberg interaction is made more isotropic, T 2 
slightly decreases. Results for 7 = 0.5 suggest that the 
spin ladders are not gpkig to be fully polarized in the 
isotropic limit 7 = l.O.EZl As discussed below, the corre- 
lations along the stripe S5 and 5*6 are much smaller than 
the previous ones and hence they show a more erratic be- 
havior, specially below the crossover temperature. Above 
this temperature, these correlations are AF, and the NN 
correlation is in general larger (in absolute value) than 
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the NNN one, although for J = 0.35 and just above T 2 
the opposite behavior is true in agreement with the anal- 
ysis made in Ref. [l(] which is valid for small J/t. A note of 
caution should be make. Since translational invariance 
in the direction perpendicular to the stripes is broken, 
the spin-spin correlations depends on the hole density on 
each site (e.g. (SfSf) = (rii)) which in turn has a smooth 
variation with temperature as shown in Fig. 0(b)). How- 
ever, the changes in the spin-spin correlations, specially 
at low temperatures, are much stronger than the varia- 
tion of the hole density, so one could safely ignore that 
dependency. 
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FIG. 3. Spin-spin correlations (defined in the text) vs. tem- 
perature on the 8x8 cluster, for (a) J/t — 0.35, e s = 1.0, 
7 = 0, (b) J/t = 0.7, e s = 0.75, 7 = 0, (c) J/t = 0.7, e s = 1.0, 
7 = 0.0 (crosses: Si, 12 x 12 cluster), and (d) J/t = 0.7, 
e s = 1.0, 7 = 0.25. (crosses: Si, 7 = 0.5). Correlations S5 
and Se have been multiplied by 5. 

Our final study concerns the other important issue that 
is the relationship between stripes and hole pairing. A 
sign of hole attraction is the presence of largest hole-hole 
correlations at smallest distances. Taking into account 
the remark made earlier about the broken translational 
invariance the hole-hole correlations C(r,-,rj) = (niUj) 
result proportional to the hole density at each site. It 



is then expected that C(Yi,Yj) along the stripe, C s (r), 
as shown in Fig. ||(a), present a smooth increase as the 
temperature is reduced, while the correlations along the 
first column next to the stripe, C±, are smoothly decreas- 
ing. However, near and below T 2 * these correlations be- 
have roughly independently of T. The same behavior can 
be observed for C(rj,rj), with rj on the center leg of a 
three-leg ladder and Yj on the column next to it, C c i, 
(Fig. ^(b)). The normalization adopted is such that: 

J2C((x,y ),(x,y)) = (N h , x ) 
y 

where yo is the y-coordinate of a reference site on column 
x, and (Nh,x) is the number of holes on that column. 
Then, C(r r) = 1. 

In Figs. 0(c) and (d) we show C(r,j, r^) at several points 
(rj,rj) after being averaged in a region AT w 0.2 above 
and below the crossover temperature. It may be noticed 
that there is no abrupt changes as T 2 is crossed since the 
averaged correlations immediately above and below it fall 
within each other error bars. The first place to look for 
hole pairs are on the stripes, where the largest hole den- 
sity is located. The results for the hole-hole correlations 
along the stripes (Fig. 0(c)), for all the parameter sets 
studied, show that C s (r) are smallest at nearest neighbor 
(NN) sites and largest at the maximum possible distance, 
although they are approximately constant beyond NN 
sites. This behavior is consistent with a metallic behav- 
ior of the stripes, as expected in the cupratesfl Results for 
J = 0.7 are virtually indistinguishable from the ones for 
J = 0.35, for the same e s = 1. In the metastable in-phase 
state, degenerate with the anti-phase state within error 
bars at the same low temperatures, the largest correla- 
tions also occur at the largest distance but with a more 
pronounced 7r-modulation, indicating a stronger coupling 
with the spin surrounding. The same short distance re- 
pulsion is obtained for C\(r) and also when belongs 
to the stripe and y^ to the first column next to it (C s i) 
(Fig. (|(c)). hi Fig. ||(d), similar results are shown fotthe 
12 x 12 cluster, x = 0.083. An alternative scenaricO as- 
sumes that hole pairs with d x 2_ y 2 symmetry are formed 
due to short-range AF correlations inside the spin do- 
mains located in between the stripes as well-known re- 
sults on small cluster calculations showEa Again, due 
to the dependence of the hole-hole correlations with the 
local hole density, the correlations between sites along 
the central leg of the spin domains are extremely small 
and hence they are almost completely masked by errors. 
However, the correlations between a site on the central 
leg of the three-leg ladder and a site on the first column 
next to it, which are almost identically zero on the 12 x 12 
cluster, acquire on the 8x8 cluster and x = 0.125, very 
similar values to those of C s \. These correlations could 
then become important as the hole density is further in- 
creased. Fig. ||(c),(d) also show exact diagonalization 
results for hole- hole correlations on L = 8 and L = 12 
t - Jz rings at quarter filling (T = T 2 * = 0.08, 0.18 re- 
spectively). They are very similar to the correlations C s 
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along the stripes except for a small shift due to the fact 
that the hole density on stripes is slightly smaller than 
x = 0.5. Besides, as shown earlier, Si shows that con- 
tiguous spin ladders are almost completely AF coupled 
in a 2D square lattice. This combined behavior suggests, 
at least for the t — Jz limit, a generalized spin-charge 
separation, or more properly, a Reparation between the 
spin background and the stripes.E3£3 
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FIG. 4. (a) Hole- hole correlations along the stripe, C s , vs. 
temperature. Lattice size, J, e 3 , and 7 are indicated on the 
plot. Open (full) circles and triangles correspond to NN (max- 
imum distance) sites. Squares and stars indicate C a i (see 
text) at the maximum distance, (b) Hole-hole correlations on 
the "spin" ladder, C c i- Open (full) symbols correspond to 
NNN (maximum distance) sites, (c) Hole-hole correlations 
for the 8x8 cluster, J = 0.35, e s — 1.0, 7 = above 
(open circles, squares and triangles) and below (full circles 
and squares, plus) TJ. C a (squares), C s i (circles), and C\ 
(triangles, plus) are shown. The last two correlations have 
been multiplied by 10. The down triangles correspond to C s 
but for the metastable in-phase state. The stars indicate ex- 
act results for an eight site t — Jz ring at quarter filling at 
T = T 2 * . (d) same as (c) but for the 12 x 12 cluster, J = 0.7, 
e s = 1.0, and 7 = 0. 



between the formation of the stripes at a charge order- 
ing temperature (more or less arbitrarily defined in our 
model) and the spin ordering process which takes place 
at a much lower temperature. This lower crossover, at 
which the spin domains become anti-phase domains and 
an incommensurate magnetic order appears, should cor- 
respond in the cuprates to the opening of the pseudogap. 
In this sense, this result that stems from a model where 
the stripes are caused by long-range Coulomb interaction, 
electron-phonon couplings, or other mechanisms that are 
described by an on-site potential, is at variance with re- 
cent results in W-hich the stripes are originated ia the 
pure t — J modelE'Ej In that and other approaches^!!, the 
stripes are the consequence of anti-phase domain forma- 
tion and both features should occur simultaneously. The 
results of the present study, including the behavior of 
spin-spin correlations along the stripes, open the possi- 
bility of experimentally discriminate the mechanism lead- 
ing to stripe formation and hence to determine to what 
extent the stripes are universal to the cuprates or depend 
on particular details of the various compounds. The hole- 
hole correlations along and near the stripes show a metal- 
lic behavior with no indications of hole attraction. The 
question arises if a hole attraction on the stripes could 
appear by taking the isotropic Heisenberg term in the 
model. One shoulcLtake into account that, as previous 
exact results show,Eil hole attraction is actually enhanced 
in the t — Jz model with respect to the fully isotropic t — J 
model. Besides, to give more support to our result, it has 
been suggested that stripes could be introduced in a uni- 
form t — J model-by taking an Ising spin interaction at 
the stripes links E3 An improvement on statistical errors 
in order to deal with larger hole densities and lower tem- 
peratures would be necessary to detect signs of hole at* 
traction inside the intervening regions between stripes.E£l 
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In summary, we have studied an anisotropic t — J 
model, close to the Ising limit, where straight site- 
centered stripes are imposed by an on-site potential re- 
flecting a mechanism which is not intrinsic to the 2D 
short-range electronic correlations of that model. The 
results of the present study suggest that we are able 
to study with QMC simulations the temperature region 
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